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The classical and quantum mechanics of isolated, nonlinear resonances in 
integrable systems with N > 2 degrees of freedom is discussed in terms of 
geometry in the space of action variables. Energy surfaces and frequencies 
are calculated and graphically presented for invariant tori inside and outside 
the resonance zone. The quantum mechanical eigenvalues, computed in the 
semiclassical WKB approximation, show a regular pattern when transformed 
into the action space of the associated symmetry reduced system: eigenval- 
ues inside the resonance zone are arranged on Af-dimensional cubic lattices, 
whereas those outside are, in general, non-periodically distributed. However, 
A'^-dimensional triclinic (skewed) lattices exist locally. Both kinds of lattices 
are joined smoothly across the classical separatrix surface. The statements are 
illustrated with the help of two and three coupled rotors. 



1. INTRODUCTION 



The surfaces of constant energy H{I) = E in the space of action variables / = (/i, . . . , In) contain 
the essential information about the dynamics of a compact integrable system, such as the fundamental 
frequencies u) = dH/dl and the foliation by invariant tori. The frequencies are given by the normals of 
the surfaces, whereas the foliation can be read off from the way an energy surface is divided into several 
patches, each representing a certain type of motion. "Simple systems" with trivial foliation like uncoupled 
harmonic oscillators have globally continuous and smooth energy surfaces. For many of the non-simple 
systems considered so far [1, 2, 3, 4, 5, 6, 7, 8, 9], the energy surfaces are nonsmooth at separatrices, but, 
nevertheless, continuous or can be made continuous by a reduction of the system's discrete symmetries. 
In [10, 11, 9] such systems are referred to as "one-component systems". 

The subject of the present work is a class of one-component systems composed of integrable approxi- 
mations of near-integrable systems of the form 



where e > is the perturbation parameter. Hi is a periodic function in the angles if, and Hq describes 
the unperturbed system which is required to be simple in order to ensure that its action-angle variables 
J, (f can be defined globally in phase space. The KAM theorem [12, 13, 14] states that a finite fraction 
of the unperturbed invariant tori survives smoothly deformed, namely those tori "sufficiently far" from 
resonance surfaces 



with relatively prime integer vectors m (an entire set mi, . . . ,mN has no common divisor). The re- 
maining fraction consists of chaotic trajectories, primary "islands" associated with the resonances of the 
unperturbed system (2) and higher-order islands. The smoothly deformed tori can be approximated by 
canonical perturbation theory. In the first order, the Hamiltonian (1) is averaged over the unperturbed 
tori, implicitly assuming that all angles ip are rapidly varying. Close to a primary island this assumption 
is not valid; the phase m • y> is almost stationary. Resonant perturbation theory [15] then suggests a 
better procedure: average over submanifolds of the unperturbed tori, m- ip — const; expand the resulting 
integrable Hamiltonian 



H{J,ip) = Ho{J)+eHi{J,v) , 



(1) 



m ■ uJo{J) — , 



(2) 



H{J,v) = Ho{J)+eV{J,m-v) 



(3) 
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at the resonance surface in J and in a Fourier series in ip, and retain only the most important contribution 

V{J, m-(p) = f{J) cos {qm ■ (p) , (4) 

where q is an integer. The regular dynamics in- and outside an isolated primary island described by 
the integrable Hamiltonian (3)- (4) is well understood; see, e.g., [16]. Hamiltonians of this type have 
been frequently used as physical models, e.g., for energy transfer in triatomic molecules; sec [17] and 
references therein. However, energy surfaces have only been presented for a special case with two degrees 
of freedom [3]. We here show the energy surfaces for a broader class of isolated-island systems with N >2 
degrees of freedom. 

Actions were the central ingredients for the old quantum mechanics before 1926. Bohr and Sommerfeld, 
among others, computed energy spectra by discretizing classical action integrals in integer multiples of 
h, Planck's constant divided by 2tt. The necessity of classical intcgrability was pointed out by Einstein, 
who formulated the quantization rules in terms of invariant tori and action variables [18]. Later, Brillouin 
derived from Schrodinger's equation that the quantization of actions is an approximation rigorously valid 
only in the (scmi-)classical limit h 0. Keller finally corrected this semiclassical approximation by 
Maslov indices in the presence of caustics [19]. According to the Einstein-Brillouin-Keller (EBK) rule, 
the quantization of an integrable system is a discretization of its action space by a A/'-dimensional cubic 
lattice with lattice constant h. Unfortunately, this plain recipe only applies to simple systems because the 
presence of separatrices destroys the applicability of the EBK rule, and, more severely, there is in general 
no one-to-one correspondence between classical action variables and quantum eigenvalues [10]. However, 
the recipe can be extended to one-component systems by introducing the action space of the associated 
symmetry reduced system [10]. It has been found in [7, 20, 11] that the classical partition of this space 
into domains of different types of motion carries over to the discretization: away from the separatrix 
surfaces there exist A^-dimcnsional cubic lattices, each related to the EBK rule for the corresponding 
type of motion; across the separatrix surfaces the lattices are smoothly connected. We here demonstrate 
that the Hamiltonian (3)- (4) gives rise to novel, less symmetric eigenvalue lattices. 

The outline of the paper is as follows. In Sec. 2, we transform the Hamiltonian (3)-(4) to a simpler 
form and expand it at the resonance surface in accordance to resonant perturbation theory. We then 
compute action variables and energy surfaces. The latter are graphically presented for a simple model, 
two and three coupled rotors. The section ends with a comparison to canonical perturbation theory. In 
Sec. 3, we firstly derive the semiclassical quantization condition, and then discuss the eigenvalue pattern 
in action space, illustrated with the help of the coupled-rotors model. Finally, we briefly draw conclusions 
in Sec. 4. 



2. ENERGY SURFACES IN ACTION SPACE 
2.1. TRANSFORMATION TO STANDARD FORM 

First of all, we transform the Hamiltonian (3) such that m — > (0, . . . , 0, 1). If only one component of m 

is nonzero then this is achieved by redefining the indices. In the general case we first redefine the indices 
such that the first two components mi , m2 are nonzero and relatively prime. Second, we introduce new 
phase space variables (-P,i?) with one of the new angles i?iv = m - (p being stationary at the resonance 
surface. For this purpose, we apply the generating function F(ip,P) = P ■ Qip of Goldstein type 2 [21]. 
The relations t? = dF/dP and J = dF/dip give 

p = (g*)-ij, ^ = Q^. (5) 

This transformation is not only canonical but also unimodular provided that the N x A/'-matrix Q has 
integer- valued components and determinant ±1 (For our purpose it is sufficient to consider matrices with 
determinant 1). This property ensures that (-P, t?) are action-angle variables of the unperturbed system. 
Remarkably, non-unimodular transformations are often used in the literature, even though the resulting 
variables are not action-angle variables in the sense of Liouville-Arnol'd [22], i.e. fixing the actions and 
varying the angles independently from to 27r does not yield a single complete cover of a torus. This has 
already been mentioned in [23] where a construction of the matrix Q has been given. We here make up 
a similar, but simpler matrix. For two-degrees-of-freedom systems we simply choose 

0.-f*-*) (6) 

y mi m2 y 
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where di , d2 are integers satisfying the diophantine equation 

det Q2 = dimi + d2m2 = 1 • (7) 

It is known from number theory (see, e.g., [24]) that for given mi, 1712, there exist a fundamental solution 
d\,d2, which can be calculated by Euclid's algorithm or chosen by hand. Prom this one gets a whole 

series of solutions di{n) = di — m2 n, d2{n) = d2 + mi n with an integer n. With the matrix Q2, the 
{N -2) X {N - 2) unit matrix 1n-2 and the (iV - 2) x 2 matrix 

W=(" ■ " ) (8) 

we construct for more than two degrees of freedom the N x N matrix 




Using this matrix in transformation (5), Hamilton's function (3) takes the new form 

H{P, i?jv) = Ho{P) + e V{P, ^n) . (10) 

The tilde ~ will be dropped henceforth. The angles . . . ,'!?Ar-i do not appear in the new Hamiltonian. 

Hence, their conjugate momenta Pi, ... , Pn-i arc constants of motion. As functions of the old momenta 
J only, they are in involution. The Hamiltonians (3) and (10) are therefore completely integrable. Action 
variables are introduced by means of 



Ii = ^f Pd^^ j = l,...N. (11) 



With the transformation (5) being unimodular, a set of fundamental paths 7^ on a given invariant torus 
is determined by "di = const, i ^ j. N — 1 action integrals simply are Ij = Pj if j < N, and it remains 
only one non-trivial integral, 

In = ^ {pn d-ffN , (12) 



2n _ 

where Pn is regarded as a function of ^?jv and the constants E,Pi, . . . , Pjv-i- 

2.2. EXPANSION AT THE RESONANCE SURFACE 

The resonance surface in P-space is given by Pjv = ^(Pi, . • . ,Pn-i), where A is implicitly defined 
through t?jv = dHo/dPN ~ 0. We expand the Hamiltonian (10) on this surface in the direction per- 
pendicular to it. For this purpose, jP/v — A\ is assumed to be of order ^/£, as usual in the analysis of 
resonances; see, e.g., [15, 16]. It is therefore sufficient to consider V as independent of Pjv when expanding 
H in powers of |Pjv — A| to quadratic terms 

H{P,^n) = Ho{Pu...,Pn-uA) 

+ '^-H'^{Pi,...,PM-i,A){PM-Af (13) 
+eV{Pi,...,PN-i,A,^N) , 

where the prime ' denotes a derivative with respect to Pn- The local approximation (13) is useful 
provided Hq ^ is fulffiled, which calls for a nonlinear dependence of H on Pjv; this is why it is called 
"nonlinear resonance" . It is illuminating to have a closer look at Hq in terms of the old coordinates 

H'^ = m.^-^m, (14) 

where the Jacobian duJo/dJ is evaluated on the resonance surface. Note that since Eq. (2) holds at 
resonance, m is a tangent to the unperturbed energy surface. Hence, Hq is the rate of frequency change 
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on the energy surface perpendicular to the resonance surface or, equahy, a measure of the; curvature of 
the energy surface in the direction of m. From this it is obvious that systems with planar unperturbed 
energy surfaces, like coupled harmonic oscillators, have to be treated separately. 
Let us consider Eq. (13) as a one-degree-of-freedom Hamiltonian 

H{Pn,^n) = ^{Pn - Af + Eo + eVi'dN) ■ (15) 

The "mass" M = 1/Hq can be positive or negative. The zero level Eq = Hq does not affect the dynamics. 
It is worth mentioning that due to the "one-dimensional vector potential" A, the Hamiltonian (15) is 
not invariant under time reversal, which is Pjv — > —Pn when the other momenta are regarded as fixed 
parameters. We perform a further simplification, usually called "resonance-centre approximation" , which 
ignores the dependence of the potential (4) upon the constants Pi, ... , Pjv-i, i-e. /(Pi, . . . , Pjv-i, ^) 
f. Absorbing |/| in the perturbation parameter £ > and using the freedom of shifting the angle 
i?jv — > 'i^AT + const, we take 



Vi-ffN) = sign(M/) cos (q^N) ■ (16) 

With En = {E—Eq) sign(M) as the energy of the A^th degree of freedom, we get from energy conservation 
F(Pjv,i?jv) =-E 

Pn = ±\/2\M\[En - e cos {q^N)\ + A . (17) 



This is the textbook planar pendulum \i q ~ e ~ 1 and ^ = 0; see, e.g., [16] and for a quantum mechanical 
treatment see [25]. Inspection of Fig. la reveals that for general parameters the phase portrait differs 
from that of a pendulum in that there is a chain of q identical islands centred at Pjv = A instead of 
just a single island centred at Pjv = 0. Let us specify briefly the invariant curves. At a fixed energy 
Em > £, there exist two invariant circles representing rotational-like motion (not necessarily related to 
physical rotations) with opposite sense of rotation A = sign(i?jv) = ±1, where iJjv = Pat — A. These 
circles are smooth deformations of the unperturbed ones. All other invariant objects, usually subsumed 
under the term "isolated resonance zone" , are created by the perturbation: a separatrix and q unstable 
(hyperbolic) equilibrium points at the critical energy E^ = S', q invariant "island circles" representing 
oscillatory motion in one of the potential wells labelled by S = 1, . . . , g at fixed En E {—e,e); q stable 
(elliptic) equilibrium points at the lowest energy En = —£■ It is to emphasize that the phase-space 
embedding of the island tori is topologically different from that of the smoothly deformed tori. 

The system exhibits 2q discrete symmetries. It is invariant under reflections with respect to '^n = (5^~ 
l)7r/g, S = 1, . . . ,g. These symmetries can be removed by implementing elastic reflections ^n —i^n 
about i?jv = and i^n = T^/q as depicted in Fig. lb. This symmetry reduction restricts the oscillations 
to one half of the first potential well, whereas it converts the rotations to oscillations within the interval 
[0,7r/g]. 

2.3. ACTIONS AND ENERGY SURFACES 

We now turn to the calculation of action variables. At fixed En > s and fixed A, the action integral (12) 
is the area in the (P/v, t?Ar)-plane between the invariant circle and the line Pat = as illustrated in Fig. la. 
We choose the closed integration path to be parametrized by -Sn increasing from to 27r. The action In 
takes on positive as well as negative values depending on the sense of rotation A = ±1 and, remarkably, 
on A. At fixed En < s, In is the area enclosed by the invariant circle. The chosen integration path 
lies inside a given potential well going from the left turning point (iJat = 0) to the right one (along 
the '0N > 0-branch) and back (along the '&n < 0-branch). In is positive and the same for all wells. 
Please, pay attention to the fact that In changes discontinuously upon traversing the separatrix. As 
a consequence, there is no unique limiting action we could assign to the separatrix and the embedded 
unstable equilibrium points. Instead, there are three different actions arising from three different energy 
limits, namely En — > e from below and En s with A = ±1 from above. This is in strong contrast to 
the continuous behaviour of the action In of the symmetry reduced system; see Fig. lb. At fixed En < e, 
the integration path goes from the left turning point to the "solid wall" at -Sn = i^/q and back after 
being refiected. With increasing energy, the left turning point wanders towards {In =0. Upon crossing 
the separatrix the smooth turning point is replaced by a refiection a,t §n =0. This does not spoil the 
continuity of In (but its smoothness) since at En > £ the action integral is the area between the two 
branches of the invariant circle (Pjv > A and Pjv < A) as illustrated in Fig. lb. Note that In does not 
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FIG. 1. a) Sketch of an isolated resonance zone with g = 2 in the (PjVi '?jv)-pla'ne; see Eq. (17). The lines i9jv = 
and i9jv = 27r are identified. Shaded regions represent action integrals. Filled circles and squares mark stable and unstable 
equilibrium points, respectively. The thick dashed line is the integration path Cq of the tunnel integral specified in Sec. 3. 
b) Symmetry reduced resonance. Dashed lines symbolize elastic refiections. 



depend on A as opposed to /jy. This is related to the remarkable fact that the symmetry reduction here 
does not only reduce phase space area by a factor, it also shifts its value by a constant. We call this a 
"non-trivial symmetry reduction" . 
The calculation of the action integrals is straightforward, and gives 



^^l^^l^^^(Vfc) if ^->^ (18) 
^^y^\[E{k)-{\-e)K{k)\ otherwise 



and 



J (T^ I^qlN + A if EM>e 

21 N otherwise, 

where JC{k) and £{k) are the complete elliptic integrals of first and second kind in the notation of [26, 27], 
with modulus fc^ — [Em /s + l)/2. Let us now come back to the iV-degrees-of-freedom system by noting 
that Ij = Ij ~ Pj a j < N (we do not care about further possible discrete symmetries). The system has 
the one-component property stemming from the continuity of In . This is mirrored by the geometrical 
fact that the two different regions of /-space, the interior of the resonance zone with En {I) < e and 
the exterior with En{I) > £, are continuously connected at the "separatrix surface" En{I) ~ e. The 
situation is more involved in /-space. The interior of the resonance zone consists of q identical parts 
labelled by £ = 1, . . . , g. The exterior is made of two parts labelled by A = ±1 which are separated by a 
gap with size proportional to the square-root of the perturbation parameter. 



A/ - 2qiNiEN = e) - V£\M\- . (20) 

TT 

A special situation occurs in the limit of vanishing perturbation strength. The resonance zone disappears, 
whereas its exterior coincides with the unperturbed action space after applying the inverse of transfor- 
mation (5). It is therefore reasonable to subject this region to the inverse transformation also for finite 
perturbation, while keeping the other region as it is. This changes the set of fundamental paths only in 
the former region, which is allowed because the other region is separated by a separatrix which prevents 
a smooth continuation of fundamental paths anyway. We thus introduce new action variables as 

J / inside the resonance zone , . 

" \ Q*/ outside. ^ ^ 

In i-space, the gap between the two parts outside the resonance zone is 



AL = \m\AI = y^e\M\\m\- . (22) 

TT 
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FIG. 2. (a) Energy surface E = 1/2 and (—1, l)-resonance surface of two free rotors in J-space; (b) transformed energy 
and (0, l)-resonance surface in P-space. 



Comparison with Eq. (14) brings to light that Ai does not depend on the length of m. Roughly speaking, 

AL is small (large) if the curvature of the energy surface in m-direction is large (small). 
The frequencies oj = dH/dl are calculated analogously as the actions giving 

oj^=s\gn{M){^'] ' (23) 



N 



and 



dEc, IT A\ UN dM ,, dA if ip ^ ^ 

Wf-lN^^ Otherwise 



(24) 



with j < N. 

2.4. EXAMPLE: COUPLED ROTORS 

Let us illustrate the previous considerations with an example of A'' coupled identical rotors described 
by the Hamiltonian 

H{J,(p) = ^J'^ +e cos (qm- (p) . (25) 

There is no need to expand this fimction; the unperturbed Hamiltonian Hq = J^/2 = J ■ J /2 is already a 
polynomial of second degree in J and the perturbation is independent of J. The space of the unperturbed 
action J is foliated by (A/^ — l)-dimensional concentric energy spheres Ho{J) = E. The resonance surfaces 
m • Wo = m • J = form a dense set of (A^ — l)-dimensional planes passing through the origin. Figure 2a 
shows such surfaces for two degrees of freedom. The energy surface provides a comprehensive picture of 
the dynamics (free motion of two particles with coordinates <pi and if2 on a circle) at a fixed energy, in 
contrast to phase portraits as in Fig. 1 which contain only information about a single degree of freedom. 
The one-piece energy surface is present in all four quadrants of action space. From this we can infer the 
existence of only a single type of motion with two rotational degrees of freedom. Each point on the energy 
surface corresponds to an invariant 2-torus (a two-dimensional torus) in phase space with the outward 
normal being the torus' fundamental frequencies. The most important information contained in the two 
frequency components is their ratio, the winding number. A rational winding number indicates that the 
periods of the rotors are rationally related. For example, a ( — 1, l)-resonance implies identical periods, 
so the particles move around the circle synchronously. In phase space, this motion is a periodic orbit. 
The time-independent phase difference <p2 — fi parametrizes a one-parameter family of such periodic 
orbits, forming a resonant 2-torus. These resonant tori are located on resonance surfaces in action space. 
For two degrees of freedom, the resonance surfaces are one-dimensional (see Fig. 2a) but we nevertheless 
refer to them as "surfaces" . Resonances are of great importance due to their sensitivity to perturbations. 
Under general perturbations, chaotic motion spreads out from resonances (and separatrices). 

For some special m with q = 1, the coupled-rotor model has a simple physical interpretation. First, if 
all numbers mi, . . . ,mjv are zero except rrij, the perturbation can be regarded as an harmonic potential 
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FIG. 3. Energy surface E = 1/2 of two coupled rotors with m = (—1, 1), q = 1, and e = 0.02 in I- (a) and £/-space 
(b). The dashed lines indicate the (0, 1)- (a) and the (—1, l)-resonance (b) of the unperturbed system. The dotted line 
serves for the construction of the quantity . 

of a spring connecting particle j with a fixed point on the circle. Note that the spring has a "negative 
spring constant" for positive e. Second, rrii = 1, mj = — 1 and all other components vanishing models 
a spring between particle i and particle j. We discuss this case for two degrees of freedom in more 
detail. In order to have a unimodular transformation (5) with m = (—1, 1), we choose (^1,^2) = (0, 1) 
leading to M = 1/2, A = Pi/2, and Eq = Pi/4:. The new angle i?i is equal to the old angle ipi and 
^2 = 'P2 — ^i describes the relative motion of the two rotors. It is worth mentioning that the most intuitive 
transformation, relative coordinate 'd2 and "centre of mass" coordinate 'di = (pi + (p2 or i}i = (tpi + V2)/2, 
is not unimodular. Yet, it is essential that the transformation (5) is unimodular, otherwise the normals 
of the transformed energy surface plotted in Fig. 2b would not give the fundamental frequencies of the 
motion on tori. Note that the (— 1, l)-resonance (^2 — 'Pi — 0) is transformed into a (0, l)-resonance 

(l?2 = 0). 

Both energy surfaces in Fig. 2 are borderless which is quite an untypical feature in the class of systems 
studied so far. Ordinary energy surfaces of two-degrees-of-freedom systems consist of different patches 
bounded by critical points. The critical points are related to isolated periodic orbits, indicating bifurca- 
tions of invariant tori. A critical point is called elliptic or hyperbolic depending on whether the periodic 
orbit is stable or unstable [8]. At an unstable orbit which is always accompanied by a separatrix the 
energy surface has a singular curvature at the critical point. It is natural that one action is zero at a 
stable orbit [3]. This can be achieved by a proper choice of fundamental paths on the invariant tori. 

The perturbed energy surface in /-space is calculated according to Eqs. (18)-(19) by fixing the energy E 
and varying the momenta Pi, ... , P^v-i- The result shown in Fig. 3a is an energy surface which consists 
of four patches and is more generic than the unperturbed one. The isolated resonance zone appears 
with two small symmetric patches with > /" = 2y/E — e which are related by time reversal. A 
point on these patches belongs to an island torus where the old angle t?2 describes oscillations, so only 
positive values of I2 are meaningful. The two hyperbolic points with (|/i|, I2) = [I^ , A/ = \/2e4/7r) mark 
unstable periodic orbits, clockwise and anti-clockwise rotating with 1^2 = ^^2 = 0, and separatrix motion, 
which is asymptotic to the embedded unstable periodic orbit. The two elliptic points with maximum 
and /2 = characterize stable periodic orbits, clockwise and anti-clockwise rotating with -82 = tt, i?2 = 0. 
Outside the resonance zone there exist two patches provided that E > e. If i? ^ e as in Fig. 3a, both 
patches together look like the unperturbed energy surface in P-space shown in Fig. 2b, apart from the gap 
of size A/ where the unperturbed surface has the (0, l)-resonance. Points on these patches correspond to 
rotational motion similar to the unperturbed motion with both old angles di and i?2 covering the entire 
interval [0, 27r). The perturbation just lifts the Constance of the velocities ??i and ■02- The patches are 
bounded by two pairs of hyperbolic points with Ii = ±7". Each pair is related to one of the unstable 
periodic orbits and separatrices discussed above. Elliptic points do not exist. 

Figure 3b displays the energy surface after transformation (21) is applied. Its rough features are 
captured by the following slight modification of the unperturbed energy surface in Fig. 2a: cut in holes 
of size AL = a/eB/tt at the (— 1, l)-resonance surface; draw a line from the intersection point of the 
resonance and energy surfaces tangential to the energy surface as pictured in Fig. 3b; its intersection 
point with the Li-axis, L'', is related to via 



L' = \m2L^ 



(26) 




FIG. 4. Winding number of two coupled rotors with E = 1/2, m = (—1,1), and q = 1. The symmetric regime 
of negative actions is omitted, a) Winding number Wo{E,Ji) of the unperturbed system. The dotted line indicate the 
resonance Wo{E,Ji) = 1. b) Wo{E,Pi). c) W{E,Ii) with e = 0.02. Inset: neighbourhood of the separatrix Ii = d) 
W{E,Li). 



add two, almost vertical pieces of height AJ at |Li | = L** = In our example with m2 = 1 the quantities 
and L'^ are equal. In general, > L'^. 

The fine structure of the perturbed energy surface is illustrated in Fig. 4 with the help of the wind- 
ing number. The unperturbed winding number Wq = ^2/'Pi is simply ±-\/2£' — J\l J\. The (—1,1)- 
resonance is characterized by Wo(-B, Ji) = 1. The same resonance in P-space is given by Wq{E^P\) = 
'^2/'^! = 0. Comparison of Figs. 4b and 4c shows that the winding number W{E,Ii) at finite e does 
not differ much from Wo{E, Pi), apart from a new piece inside the resonance zone, which takes on small 
values of order ^/s. The behaviour at the separatrix Ii = P can be seen more clearly in the magnification. 
Exactly at the separatrix, W{E,Ii) logarithmically approaches zero, or, taking a more common point of 
view, l/W{E,Ii) diverges logarithmically. This means, on the one hand, that the associated unstable 
motion can be regarded as resonant. On the other hand, it means that even though there is a gap in 
the energy surface, there is no such gap in the spectrum of the winding number. It can be inferred from 
Fig. 4d that this is also true for the transformed winding number W{E,Li). It approaches the value 1 
(in general —m2/m\) at the separatrix. The derivative of the winding number is large in the vicinity of 
the separatrix. As a consequence, there is an accumulation of low-order resonances {W is a fraction of 
two integers with small denominator) near the separatrix. 

Figure 5 shows how the L-space is foliated by energy surfaces. The evident discontinuity differs strongly 
from the continuity of J-space, cf. Figs. 5 and 6. Let us try to get more familiar with the symmetry 
reduction on the basis of the model of two rotors coupled by a spring, m = (—1, 1) and q = I- The 
symmetry-reducing reflections about = and = tt are related to two symmetry transformations, 
^ —'&2 and i?2 — TT — > TT — i?2 • The first one reads in the old angles: <^2 — </5i ^ '^i — '^2- This is an 
interchange of particle 1 and 2, so the reflection at i?2 = can be viewed as an elastic reflection between 
both particles. The interpretation of the second symmetry transformation, (^2 — {'■Pi+t^) (</'i +7'') — V2, 
is more involved: shift particle 1 by tt on the circle; interchange both particles and flnally shift particle 1 
by — TT. The reflection at '§2 = 'k can be seen as an elastic reflection of particle 2 with particle 1 virtually 
displaced by tt. 

Let us take a brief look at an example with m2, 9 > 1, namely m = (—2, 3) and q = 2 with (di, c?2) = 
(1, 1). We observe from Fig. 7 that the energy patches are well separated as predicted by Eq. (26). A point 
within the resonance zone represents two 2-tori, each a combination of a i^i-rotation and a ■!?2-oscillation 
in one of the two potential wells. 

The description of the two-degrecs-of-frecdom dynamics extends to three (and more) degrees of freedom 
in a natural way. Wc demonstrate this with three coupled rotors, using m = (1,1,1), q ~ 1 and 
{di,d2) = (0,1). Figure 8 shows that the mipcrturbed energy surface is made of a single piece being 
present in all octants of action space, so only one type of motion exists with three rotational degrees of 
freedom. Again, the surface is somewhat special in that it is borderless. Typical energy surfaces of three- 
degrees-of-freedom systems are composed of several patches bounded by critical edges. Interior points 
of a given patch represent 3-tori. Elliptic edges correspond to stable isolated 2-tori, whereas hyperbolic 
edges correspond to unstable isolated 2-tori and separatrices. Corner points represent isolated periodic 
orbits and separatrices. Although isolated low-dimensional tori do not exist in the free-rotor model, there 
are families of low-dimensional tori on resonance surfaces. Figure 8 reveals that these surfaces intersect 
a spherical energy surface in great circle meridians. All intersection lines together form a dense set, the 



RESONANCE ZONES IN ACTION SPACE 



9 






FIG. 7. Energy surface _E = 1/2 of two coupled rotors with m = (—2, 3), g = 2, and e = 0.02. The dashed Une marks 
the (— 2, 3)-resonance of the unperturbed system. The dotted line serves for the construction of the quantity . 
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FIG. 8. Energy surface -E = 1/2 of three free rotors in J-space. Intersection lines with (mi, m2, m3)-resonance surfaces, 
\mj I < 1, are shown. 

so-called Arnold web. Along intersection lines we find resonant 3-tori, one-parameter families of 2-tori. 
A torus at an intersection point of two such lines is completely resonant, i.e. it is foliated by periodic 
orbits. The resonance surfaces are shown in Fig. 8 up to order one, i.e. \mj\ < 1. Note that the order 
so defined is not invariant under unimodular transformations. The resonances on the energy sphere in 
Fig. 8 and the transformed surface in Fig. 9 are therefore not always related via transformation (5). 

Figures 10 and 11 illustrate how the perturbation modifies the energy surfaces in 3D action space. Let 
us concentrate on the X-rcpresentation in the latter figure. The resonance zone cylindrically surrounds the 
other parts of the energy surface. Its elliptic edge in the plane L3 = presents isolated stable 2-tori. Such 
a torus is essentially a direct product of a circle and the stable periodic orbit discussed for two degrees 
of freedom. The other edge is hyperbolic and (0, 0, l)-resonant. Its unstable 2-tori and separatrices are 
again direct products of a circle with the corresponding two-degrees-of- freedom object. The resonance 
zone has no corner points and, correspondingly, no isolated periodic orbits. The energy surface outside 
the resonance zone is torn open along the (1, 1, l)-resonance surface of the unperturbed system; cf. 
Figs. 8 and 11. Apart from this gap, the shape of the surface is again similar to that of the unperturbed 
surface. The hardly visible deformation close to the gap is uncovered by the graphical representation of 
resonances. The two hyperbolic edges are special (1, 1, l)-resonances which are tangentially approached 
by other resonances. This is analogous to the earlier- mentioned accumulation of resonances in the vicinity 
of separatrices in the case of two degrees of freedom. 



RESONANCE ZONES IN ACTION SPACE 



11 




FIG. 9. Energy surface E = 1/2 of three free rotors in f'-space. 




FIG. 10. Energy surface E = 1/2 of three coupled rotors with m = (1, 1, 1), g = 1, and e = 0.02 in /-space. The 
unmarked resonances can be read off from the unperturbed energy surface in Fig. 9. 
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2.5. COMPARISON WITH CANONICAL PERTURBATION THEORY 

Having calculated the energy surfaces of coupled rotors exactly, wc now compare them to approximated 
surfaces. Canonical perturbation theory applies near-identity transformations parametrized by e to the 
action-angle variables such that the new Hamiltonian is independent of the angles and therefore integrable; 
see, e.g., [16]. After eighth-order canonical perturbation theory the Hamiltonian (25) transforms to 

H^ML) = -L^ + e^-e + + e*^— + e^^^^^ , (27) 

dppv ; 2 4^ 64^ 128^ 16384^ ' ^ ' 

with the new actions L and the abbreviation ^ = \m\/{m ■ L). It turns out that _ffapp(i) agrees with 
the Taylor scries of H{L) for E'jv > £ which is obtained from expanding Eqs. (18)-(19) in powers of 
£, transforming according to Eq. (21), and solving (to eighth order) for the energy. Hence, canonical 
perturbation theory approximates the energy surfaces outside the resonance zone, cf. Figs. 3b and 12. 
But notice that Hg,pp{L) diverges on the resonance surface of the unperturbed system, m ■ L = 0. This 
is responsible for the fact that the approximate energy surfaces do not have hyperbolic boundaries but 
instead additional segments close to the resonance surface without any physical interpretation! 

The energy surfaces within the resonance zone cannot be approximated in this way. The reason is that 
the near-identity transformations cannot cope with the different topology of the phase-space embedding 
of island tori; see, e.g., [16]. We arrive at the same conclusion by observing that for < £ the action 
/jv in Eqs. (18)-(19) cannot be expanded in a Taylor series in powers of e. 



3. DISCRETIZATION OF ACTION SPACE 
3.1. SEMICLASSICAL QUANTIZATION CONDITION 

It is instructive to begin with the EBK rule of the unperturbed system, 

J={k + (3/4)h, (28) 

with quantum numbers k and Maslov indices (3. The latter are classical indices characterizing the motion 
on a given invariant torus: (3j = if the jth degree of freedom is of rotational type and (3j 2 if the 
jth degree of freedom is of oscillatory type [19, 28]. For rotational motion, Jj and kj assume positive as 
well as negative values (example: angular momentum in a system with rotational symmetry), whereas for 
oscillatory motion, both numbers are usually restricted to non-negative values (example: one-dimensional 
harmonic oscillator). A Hamilton operator Hq with semiclassical eigenfunctions 



V'fc(¥',e = 0)cxexp[i(fe + /3/4)-<p] (29) 

can bo obtained from the unperturbed Hamilton's function Hq with the usual operator replacement 
Jj —ihd/d(pj, j — 1, . . . ,N and = —1. The semiclassical eigenvalues of Hq are E = Hq{J) with J 
from Eq. (28). Note that the functions (29) are in general not 27r-periodic in the angles ^Pj. This stems 
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from the singularities in the transformation from the original Cartesian coordinates x and momenta p to 
action-angle variables. Transforming the functions (29) according to Eq. (5) results in 



e = 0) a exp[i(n + a/4) • i?] , 



(30) 



with new quantum numbers n = {Q~^Yk and Maslov indices a = {Q~^Yf3. We see here again that it is 
important to employ a unimodular transformation; it guarantees that nj and aj are integer-valued and 
that consecutive values of rij differ by 1. 

The EBK rule cannot be applied to the perturbed system (3)- (4) because the actions and Maslov 
indices are not globally defined for the entire phase space; their definition is different for the interior and 
the exterior of the isolated resonance. We need a quantization condition which is uniformly valid for 
the entire phase space. Our derivation of this uniform quantization condition is similar to the derivation 
in [29] but differs in three respects: an arbitrary number of degrees of freedom is allowed; the potential (4) 
has to be independent of Pi, ... , Pjv-i; action-angle variables are used at every stage. Let us start from 
the ansatz 



ipni'0) oc exp 



N-l 



iJj{-dN,nN) , 



with the eigenfunction of the A^th degree of freedom 

-ipi'&NjTiN) = exp[i(ajv/4)z?jv] ri(t?Ar,njv) . 



(31) 



(32) 



The function fli^N ,71^) is to be 27r-periodic in -Sn in order to ensure Q{'dj^,nN) — > exp(mArt!?jv) and 
'4'n{'&) V'n(t?,e = 0) as £ ^ 0. Clearly, N — 1 degrees of freedom in this variables are quantized a la 
EBK 



l3=Pj = {nj+aj/4)h, j = l. 



,N-1 



(33) 



In order to find a quantization condition for the remaining degree of freedom we insert the conditions (33) 
and P/v —ihd/d'&N in Hamilton's function (15) leading to a Hamilton operator with eigenvalue equa- 
tion 



1 

2M 



-ih 



(I'd 



N 



A] +Eo-E + eV{'dN) 



■tp{'dN,nN) = 



With the Bloch-wave ansatz 

tlj{'&N, un) = ex.p{iA'&N/n) *('i?jv, njv) 
a Schrodinger equation without "vector potential" A is obtained 



'2Md'd% 



+ Eo-E + sV{'&n) 



*(i?iv,njv) = . 



The wave functions ^'(^jv,njv) fulfill "twisted boundary conditions" 

^{i^N + 2-77, njv) = exp(27ricr) *(^?jv, tin) 

with the real quantity 

£7 = -A(7i,...,7jv-i)/?i + ajv/4 , 
and d^/d'&N continuous. With the potential (16) (|/| is again absorbed in e) we finally get 

d^ 



'2\M\ d-ffl 



Ej<f -|-£cos(g^?jv) 



N 



*(^jv,njv) = . 



(34) 



(35) 



(36) 



(37) 



(38) 



(39) 
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According to [30] , the semiclassical solutions En of this eigenvalue problem with boundary conditions (37) 
are given by 

cos\2Tr(l — a)/q] , ^ 

°"'^-'''° ^i + eUe/Ir 

Let us specify the constituents of this formula. The first one is the phase integral 

(piEN) = ^ [ PncIi^n . (41) 

At fixed En < e, the integration path connects the turning points inside a potential well. At En > s, 
the path goes from i?jv = to z?iv = 27r/g. It is an easy exercise to show that the phase integral is related 
to the action integrals (18) by means of = 27r/]v/?i- The path lies entirely in the classically allowed 
region of phase space, in contrast to the integration path of the tunnel integral 

Q{En) = -i I Pn ddN , (42) 

which lies in classically forbidden regions. At fixed En < s, the path Cq links turning points from 
neighbouring potential wells through the barrier as shown in Fig. la. In this case the tunnel integral 
of the barrier is positive. At En > £, the path connects complex turning points, which are complex 
conjugates of each other, yielding a negative tunnel integral. Again we show directly the outcome of the 
calculation: 

Q ^ I '-fVmk [^(^/F^/fc) - ^^/W^/k)] for En > e ^^^^ 
\ ^^J\M\ - k'^K{^/l^)\ otherwise, 

where k is the same modulus as for the action variable (18)-(19). p is the "quantum correction function" 



P{En) = -^1o§ 



eirfi \ ' 



(44) 



with e = exp(l) and /3 = 1.78107. A detailed derivation of the formula (40) based on a WKB ansatz can 
be found in [31]. Note that even though the WKB method is only valid in the semiclassical limit % 0, 
or to put it another way, for highly excited states, it gives very often accurate results even for the ground 
state. 

For solving the quantization condition (40] it is convenient to use a combination of bisection and 
Newton's method. Firstly, the actions /i, . . . , In-\ are determined from the given numbers ni, . . . , un-x 
by virtue of the rules (33). Secondly, the minimum energy eo = Ea{Ii, . . . , In-i) — sign(M)£ with 
0(eo) = is computed. Thirdly, the running index j is initialized to 0. Fourthly, we look for an energy 
Bj+i such that the interval (cj, e^+i) includes exactly one eigenvalue. To do so, we note that the left hand 
side of condition (40) is essentially a cosine of (j) {p can be ignored for the following arguments) and that 
the right hand side and are monotonic functions of the energy. Hence, a suitable energy e^+i is given 
implicitly by the relation 0(ej+i) = jtt which is solved numerically with Newton's method taking ej as 
starting value. Having determined the interval, the bisection method is employed to find the enclosed 
eigenvalue. Finally, j is increased by one and the last steps are repeated until the desired number of 
eigenvalues is found. 

3.2. LATTICE STRUCTURE 

An analogy between action-space discretization and crystal lattices can be drawn by rewriting the EBK 
rule for the unperturbed system (28) as 

N 

J = ao + ^ kitti , (45) 

i=l 

with basis vector Oq = /3?i/4 and primitive lattice vectors a^j = 6ijTi, where 5ij is the Kronecker symbol. 
Equation (45) defines a lattice in J-space, the primitive elementary cell of which is an ^''-dimensional 
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cube with side length h. What kind of lattice does the quantization c;ondition (40) imply? Ac:cording 
to [20], we formulate the quantization condition in terms of the action variables of the symmetry reduced 
system I. This is trivial for the phase integral (/> = 27r/Ar/?i, but the relation = En{In) has to 
be computed numerically. Note that the relation is unique and continuous due to the one-component 
property; see also [10]. The tunnel integral and the quantum correction function become functions of 
In via En = En{In)- We observe that the quantization condition (40) is a function of the actions I 
alone. It is essential to realize that we cannot replace the actions of the symmetry reduced system / by 
the actions of the full system I using Eq. (19) since the classical index A is not provided by quantum 
mechanics. Quantum mechanically, we cannot distinguish between a classical torus with A = +1 and 
its symmetric partner with A = — 1. This is the reason for the importance of the action space of the 
associated symmetry reduced system. 

The quantization condition (40) has two solutions with the same quantum number I in each /jv-interval 
of width h. Wc therefore add two further quantum numbers, 11 = ±1 and njv 0, 1, 2, . . .. The former 
distinguishes both solutions and the latter label the intervals. We then verify that the quantization 
condition can be cast into the form of an EBK-like rule 

In = {fiN + aN/4:)n (46) 

by inserting this rule into Eq. (40). The quantum number tIn cancels due to the 27r-periodicity of the 
cosine. Inverting the cosine on the correct branch, which is determined by H, gives 2q "Maslov phase 
functions" 



CKjv = -arg cos[27r(l - a)/q] + exp(2@/h) + sin'^[27r(; - a)/q] ] +-p , (47) 

TT \ ' J TT 

where "arg" extracts the polar angle G [0, 2tt) of a complex number. In contrast to the type of classical 
motion and its characterizing Maslov index a^v which both change discontinuously at the separatrix, 
the Maslov phase function varies smoothly across the separatrix. Hence, Eqs. (33) and (46)- (47) do not 
define a periodic lattice in the entire action space, instead they define a "WKB lattice" in the terminology 
of [11]. 

Let us consider first the WKB lattice deep within the resonance zone. En <^ £■ From ^ h and p « 

follows an EBK rule with constant djv = 2 — 11. The /-independence of djv manifests itself in a q-iold 
quasi-degeneracy of the eigenvalues In', I = ■ ■ ■ ,q labels those eigenvalues. Translated into action-space 
geometry, the EBK rules for 7i, . . . ,7jv give two sets (11 = ±1) of g' identical lattices, each set having 
unique basis and lattice vectors. It is possible to combine both sets to a single lattice with a non-primitive 
elementary cell, a so-called "quantum cell" invented in [11]. Here, the quantum cell is a Af-dimensional 
cube containing 2q quantum states. Going back to the full system, one finds from Eqs. (19) and (46) 

Jjv = 2/Ar = {uN + l/2)n (48) 

with njv = 2njv ifn = -|-1 and un = 2njv -|- 1 if H = — 1. We see therefrom that H is the parity with 
respect to the potential's symmetry line i?jv = T^/q- The Maslov index 2 in the EBK rule (48) is in 
agreement with the oscillatory character of the motion of the A''th degree of freedom in this phase space 
regime. 

The other extreme case. En e, well outside the resonance zone, coincides with the unperturbed 
limit; we have 6 <C — ?i and p « giving 

- 4 - H 

«5;jv = n-(Z - cr) with n=^ — , . rr^ ■ (49) 

q^ ' sign{sin[27r(/-(7)/g]} ^ ^ 

Eigenvalues In with different n are twofold degenerated in the non-generic situation of 2(7 being an 
integer. The Maslov phase ajv depends on the actions Ii, . . . ,/jv-i via a{A) reflecting the non-trivial 
symmetry reduction. Clearly, the quantization conditions (33) and (46) with Eq. (49) give, in general, a 
non-periodic eigenvalue distribution in the action space of the symmetry reduced system. Nevertheless, 
we will see in the following that periodic lattices in classically small regions of action space still exist. 
Each such region, even though classically small, contains many eigenvalues in the semiclassical limit. 
Here, a is approximately a linear function of the actions 

N~l 
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Reformulating the quantization conditions (33) and (46) with the help of Eq. (50) gives an equation 
which compares to Eq. (45) (with I and n instead of J and k) describing 2q lattices with basis vectors 
aoi = aih/A and lattice vectors aij = SijU — djNllbi'h/q ifi < N, else oqat = 'n{l — bo)h/q and a^j = S^jJi- 
For fixed n and I, the lattice vectors have, in general, different lengths and arc not orthogonal; the lattice 
is iV-dimensional triclinic. Note that, in general, a single quantum cell cannot be defined due to non- 
matching lattice vectors. The reader should realize that the unperturbed eigenvalues of both the full 
and the symmetry reduced system lie on a simple A^-dimensional cubic lattice in the respective "correct" 
action space: for the symmetry reduced system, the quantum analog of the classical elastic reflections are 
Dirichlet boundary conditions, i.e. vanishing wave function, on the lines iJjv = and iJat = t^/q leading 
to an EBK rule with Maslov index 4 (a hard wall instead of a smooth turning point increases the Maslov 
index by 1). For the full system, we get an EBK rule from Eqs. (46) and (49) using Eq. (38) 

In = MIn + A= [hn + aN/'i)n (51) 

with the identifications n = —A and rijv = —AqfiN — I- In the unperturbed case, the quantum index 
n therefore equals the classical index —A which characterizes the parts outside the resonance zone. It 
is thus possible to relate each quantum wave function '5 uniquely to the classical A = +l-region or the 
A = — 1-region. At finite En > e, there is no unique relation between quantum states and these classical 
regions as already discussed. When representing the eigenvalues of the exterior of the resonance zone in 
the action space of the full system (J- or L-spacc), we project into the A = +l-rcgion (we could also 
choose the A = —1-region). In this "reduced action space" the eigenvalues lie on A^-dimensional cubic 
lattices; see Eqs. (33) and (51). 

3.3. EXAMPLE: COUPLED ROTORS 

We return to the example of two coupled rotors (25) with m = (—1, 1). The Maslov indices of free 
rotors, Pj, vanish, such as the transformed ones, aj, do. Combining with A = Pi/2 we get a = —Ii/{2Ti) = 
— ni/2. This linear dependence makes for globally periodic lattices outside the resonance zone in /-space. 
These lattices with basis vectors ao = {0,fUfi/q) and lattice vectors ai = {fi,flh/{2q)), a2 = (0,?i) are 
skewed. Because of the integer- valuedness of 2a, the lattices with II = -|-1 and the ones with n = — 1 are 
congruent, or in other words, the eigenvalues are twofold degenerated. 

Before studying the model with g = 1 in detail, it is interesting to mention that this model has already 
been treated in Born's 1925 book Vorlesungen iiber Atommechanik [32] with the old Bohr-Sommerfeld 
quantization rules (EBK without Maslov indices). This is a too crude approximation as it will become 
apparent in the following. Consider the right hand side of the quantization condition (40). It is positive 
or negative, depending on whether rii is even or odd. Figure 13 illustrates the former situation (for 
odd ni the non-periodic curve is reflected at the zero line). The action I2 is discretized according to 
Eq. (48) in the small- 72 domain (deep within the resonance zone) and according to Eqs. (46) and (49) 
with twofold degeneracy in the large- 72 domain (well outside the resonance zone). The transition between 
these different kinds of discretization happens in a narrow region around the separatrix with width of 
order h. 

Solving quantization condition (40) for all values of ni gives the complete set of eigenvalues. Figure 14 
shows their arrangement in the space of the constants of motion E and Ii = Ji + J2- The twofold 
quasi-degeneracy at large energies is related to the exact degeneracy of the unperturbed eigenvalues 
E = (kl + fc|)?i^/2 and 7i = (fci + k2)h. The eigenvalue pattern looks rather regular due to the fact 
that one constant of motion is an action variable. However, the underlying regular structures are more 
transparent in the /-representation in Fig. 15. The eigenvalues are located on the WKB lattice, which 
reduces to periodic lattices far away from the separatrix surface. Below the separatrix, i.e. inside the 
resonance zone, the eigenvalues lie on two quadratic lattices. It is here trivial to see how a larger square 
elementary cell (the quantum cell) could combine both lattices to a single one. Above the separatrix, i.e. 
outside the resonance zone, twofold degenerated eigenvalues lie on two different but congruent skewed 
lattices, the elementary cells of which are illustrated in Fig. 15. The integcr-valuedness of 2(7 ensures 
here the existence of a quantum cell, a body-centred square in the terminology of crystallography. As the 
separatrix surface is crossed from above, a smooth, degeneracy-lifting transition to the quadratic lattices 
takes place. The lattice outside the resonance zone is even simpler when projected into the A = +l-rcgion 
of i-space as displayed in Fig. 16. Away from the separatrix surface, twofold-degenerate eigenvalues are 
arranged on quadratic lattices according to the EBK rule of the unperturbed system (33) and (51). The 
price to pay for recovering a simple lattice for a subset of eigenvalues is the loss of the coherent picture 
of the action-space discretization as shown in Fig. 15. 
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FIG. 13. Graphical solution of the quantization condition (40) for two coupled rotors with m = (—1, 1), q = 1, £ = 0.02, 

and ni even. Eacii side of Eq. (40) is drawn as a function of I2 = (bTi/(2n) in units oih = 0.02; intersection points indicate 
quantized values of Iz- The dashed Une marks the separatrix I2 = A//2. 
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energy E and action I\ = J\ + J2 (in units of ?i = 0.02). The symmetric region of negative 7i is omitted. 
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FIG. 15. Semiclassical eigenvalues of two coupled rotors with m = (—1, 1), q = 1 and e = 0.02 in /-space in units 
of h = 0.02. The symmetric region of negative action Ii is omitted. The solid line is the energy surface E = 0.15 and 
the dashed line is the separatrix surface. Filled regions represent elementary cells belonging to n = +1 (left) and n = — 1 
(right). 
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FIG. 16. Semiclassical eigenvalues of two coupled rotors with m = (—1, 1), q = 1 and e = 0.02 in reduced i-space in 
units of ?i = 0.02. Eigenvalues belonging to the resonance zone are not shown. The solid line is the energy surface E = 0.15 
and the dashed line is the separatrix surface. 
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FIG. 17. Graphical evaluation of Eq. (40) for two coupled rotors with m = (—1, 1), g = 2, e = 0.02, h = 0.02, and 
ni/2 even. The dashed line marks the sepaxatrix. 

Let us now discuss the case q — 2. The quantization condition (40) simpUfics to an EBK rule if its right 
hand side vanishes. This is the case if ni is odd since a = —ni/2. In the more comphcated case of rii 
even, we have to distinguish between ni/2 even and ni/2 odd. Figure 17 illustrates the former situation. 
Note that the sign of the right hand side of the quantization condition is determined by I. (For odd ni/2 
the solutions to / = 1 and I = 2 are interchanged.) We see here in addition to the twofold degeneracy 
for large I2 also a twofold degeneracy for small h- Figure 18 shows that above the separatrix surface 
the lattice in action space looks similar to g = 1 in Fig. 15. However, the situation now is actually a 
bit more involved: two kinds of elementary cells labelled by S = ±1 are shifted by h/2 in /2-direction 
(Z = 1, 2); we have four lattices instead of two. Again, a quantum cell could be defined. A more striking 
difference between the case g = 1 in Fig. 15 and g' = 2 in Fig. 18 is that in the latter case there is a 
twofold degeneracy below the separatrix surface. 

Finally, we illustrate a more generic case with irrational a using two coupled non-identical rotors as 
example 

Ho = lj! + lji, (52) 

with 7 being the reciprocal of the golden mean, (V5 - l)/2. We again take m = (—1, 1) and q = 1 
giving £7 = — ni/(l +7). Figure 19 confirms that the quadratic lattices below the separatrix surface are 
as in the case of integer- valued 2a pictured in Fig. 15. But above the separatrix surface, there are two 
non-congruent skewed lattices with basis vectors — (0,n7i) and lattice vectors ai = (?i,n?i/(l + 7)), 
a2 = (0, ?i). A quantum cell does not exist in this case, which is expressed by the seeming irregularity 
of the overlap of the lattices in Fig. 19. Note that the eigenvalues in Fig. 19 are distinguished by their 
quantum number 11 and not by the lattice index fl. 

4. CONCLUSIONS 

The classical and quantum mechanics of isolated, nonlinear resonances has been presented in action 
space. The energy surfaces were found to be typically composed of two different kinds of patches corre- 
sponding to the inner and the outer part of the resonance zone. The graphical representation of these 
smfaces for a model of coupled rotors proved to be a very concise description of the integrable dynamics. 
Moreover, it was demonstrated that energy surfaces outside the resonance zone can be approximated 
by canonical perturbation theory, but have to be handled with care, since non-physical segments are 
produced close to the resonance zone. 
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FIG. 18. Semiclassical eigenvalues of two coupled rotors with m = (—1, 1), q = 2, e = 0.02, and h = 0.02 in /-space. 
The solid line is the energy surface E = 0.1b and the dashed line is the separatrix surface. 




FIG. 19. Semiclassical eigenvalues of two coupled non-identical rotors with irrational cj, m = (—1, Q = ^ = 0.02, 
and h = 0.02 in /-space. The dashed line marks the separatrix surface. Filled regions represent elementary cells belonging 
to n = -hi (left) and n = -1 (right). 
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Exploiting the one-component property of the systc;ni, we have investigated the distribution of the 
quantum mechanical eigenvalues in the action space of the symmetry reduced system: the eigenvalues 
within the resonance zone are located on A''-dimensional cubic lattices, whereas the other eigenvalues are 
located on locally A'^-dimensional triclinic lattices reflecting the non-trivial symmetry reduction. Both 
kinds of lattices and the smooth transition between them were described by a uniform semiclassical 
quantization procedure and graphically illustrated with the help of the coupled-rotor model. It was 
found that simple cubic lattices can be recovered only separately for each of the different types of classical 
motion in a properly reduced action space of the full system. 

The presented discussion deals with nonlinear resonances. The important linear case of coupled har- 
monic oscillators demands special considerations in a future publication. 
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